Resonant tunneling of interacting electrons in a one-dimensional wire 
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We consider the conductance of a one-dimensional wire interrupted by a double-barrier structure 
allowing for a resonant level. Using the electron-electron interaction strength as a small parameter, 
we are able to build a non-perturbative analytical theory of the conductance valid in a broad region 
of temperatures and for a variety of the barrier parameters. We find that the conductance may 
have a non-monotonic crossover dependence on temperature, specific for a resonant tunneling in an 
interacting electron system. 
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The phenomenon of resonant tunneling is well-known 
in the context of electron transport physics [1]. The hy- 
bridization of a discrete state localized in the barrier with 
the extended states outside the barrier may strongly en- 
hance the transmission coefficient for electrons incident 
on the barrier with energy matching the energy of the 
localized state. For a single electron, the transmission 
coefficient at energies close to the resonance is given by 
the Breit-Wigner formula [1]. However, if the barrier 
carrying the resonant level separates conductors which 
in equilibrium have a finite density of mobile electrons, 
the problem of resonant tunneling becomes more complex 
due to the electron-electron interaction. Manifestation of 
resonant tunneling in the conductance of a solid-state de- 
vice is inevitably sensitive to this interaction. 

Some of the effects of electron-electron interaction do 
not depend on the dimensionality d of the conductors- 
leads separated by the barrier. For instance, the on-site 
repulsion together with the hybridization of the localized 
state with the states of continua lead to the Kondo effect 
in the transmission across the barrier [1,2] at any d. The 
Fermi-edge singularity also strongly affects the resonant 
tunneling [3,4] in any dimension. The electron-electron 
interaction within the leads, however, does not have a 
strong effect if d > 1, and if the leads are not disordered. 
In contrast, resonant tunneling across a barrier interrupt- 
ing a one-dimensional (Id) wire is modified drastically by 
the interaction within the wire. The importance of such 
a setting is emphasized by the recent transport experi- 
ments with nanotubes and nanowires containing a quan- 
tum dot [5,6]. 

The electron-electron interaction enhances the 
backscattering off the barrier for electrons with energy 
close to the Fermi level [8]. We find that if the dis- 
crete level is not perfectly aligned with the Fermi level 
in the leads, or the barrier structure has even slight 
geometrical asymmetry, then the low-temperature lin- 
ear conductance decreases to zero with the temperature, 
G{T 0) 0. At sufficiently high electron energies 
(measured from the Fermi level) the enhancement of 
backscattering due to the interaction should get weaker. 



and the conventional behavior of resonant tunneling 
prescribed by Breit-Wigner formula may be restored. 
Consequently, the conductance G(T) may increase with 
the temperature being lowered. How to match these 
two opposite tendencies? We answer this question below 
by finding the proper non-monotonic crossover function 
G{T) for an arbitrary asymmetry of the barrier and ar- 
bitrary position of the resonant level with respect to the 
Fermi level, in the limit of weak interaction. The asymp- 
totes of G(T) agree with those found in Ref. [7,8] in the 
context of the Luttinger liquid theory. The universality 
of latter results was recently doubted in the theoretical 
part of [6]. We find no ground for such doubts. 

We will use an analogue of the renormalization method 
developed in [9]. Within this method, the complicated 
picture of many-electron transport is considered within 
the traditional Landauer-Biittiker elastic scattering for- 
malism. The role of the interaction is to renormalize 
the elastic scattering amplitudes. The renormalization 
brings about an extra energy dependence of these ampli- 
tudes. It was shown in [9] that in the limit of weak inter- 
action the most divergent terms in perturbation theory 
indeed correspond to the purely elastic processes, thus 
justifying the method. The advantage of the method 
is that it allows one to investigate quantitatively the 
crossover between the limits of weak tunneling and full 
transmission across the barrier. 

The original formulation of the method [9] disregarded 
the energy dependence of scattering amplitudes in the 
absence of interaction. While valid for a generic case of 
a single scatterer, this assumption obviously fails to de- 
scribe the resonant tunneling. To circumvent this, we ex- 
tend the method to arbitrary energy dependence of scat- 
tering amplitudes. First step is to derive the first-order 
interaction correction to scattering amplitudes. This can 
be readily done along the lines of Ref. [9] . The correction 
to transmission amplitude reads 
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Here the ^/^(ij) arc the reflection amplitudes for electrons 
incoming from the left (right), and the coefficients aL(i{) 
represent the interaction within the left(right) part of 
the Id wire; energies e and e' are measiircd from the 
Fermi level. Transmission and reflection amplitudes rL,R 
satisfy the unitarity relation: trI* = —r'^t. The coefE- 
cients a are related to the Fourier components V{k) of 
the corresponding electron-electron interaction potential 
by a = {V{0) - V{2kF))/2-KVF). In the limit of weak 
interaction, these coefficients determine the exponents in 
the edge density of states [8] for each part of the channel, 
v{e) oc e". 

The integration over e' in the first-order correc- 
tion Eq. (1), in general, yields a logarithmic divergence 
at e ^ 0. This indicates that the perturbation series 
in the interaction potential can be re-summed with the 
renormalization method. To account for the most diver- 
gent term in each order of the perturbation theory in a, 
we proceed with the renormalization in a usual way [10]. 
On each step of the renormalization, we concentrate on 
the electron states in a narrow energy strip around —E, 
with E > being the running cut-off. We evaluate the 
interaction correction due to the electrons in these states 
to the scattering amplitudes at energies e close to Fermi 
level, |e| < E. These amplitudes are thus functions of 
both e and E. We correct those amplitudes according 
to Eq. 1, reduce the running cut-off by the width of the 
energy strip, and repeat the procedure. This yields the 
following renormalization equation; 
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provided that |e| < E. We abbreviate here r{e) = r(e, |e|) 
(and similar for t) indicating that the renormalization 
of scattering amplitudes stops when the running cut-off 
approaches The initial conditions for this differen- 
tial equation are set at upper cut-off energy A. If the 
e-dependence of the transmission amplitude in the ab- 
sence of interaction, t(e. A), can be disregarded, then all 
the energy dependence of renormalized amplitudes comes 
about as a result of the renormalization procedure. The 
corresponding simplification of Eq. (2) then reads 
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and contains the transmission probabilities only. This 
coincides with the results of Ref. [9]. However, the above 
simplification is not possible in the more general case 
we consider here. One can not even deal with a single 
equation: the equation (2) shall be supplemented with a 
similar equation for one of the reflection amplitudes. 
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To describe resonant tunneling, we consider a com- 
pound scatterer made of two tunnel barriers with tun- 
nel amplitudes ti^2 ^ 1 separated by a distance wvp/S. 
This gives rise to a system of equidistant transmission 
resonances separated by energy d. We assume that one 
of the resonances is anomalously close to Fermi energy 
and concentrate on this one disregarding the others. The 
scattering amplitudes in the absence of interaction are 
then given by common Breit-Wigner relations: 



t(e,A) = 
r-z,(e,A) 
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where Tl,r = |ti,2p^/27r are the level widths with re- 
spect to the electron decay into the left (right) lead and 
A is the energy shift of the resonance with respect to the 
Fermi Level; we assume here A <i; 5. We disregard pos- 
sible energy dependence of ti_2 that could be relevant at 
higher energies, which allows us to take the upper cut-off 
A to be of the order of S. The corresponding transmission 
probability before the renormalizations. 
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is the usual Lorentzian function of energy. The interac- 
tion corrections to A and Tl jj which come from bigger 
energy scales, 6 < E < Ep, are assumed to be included 
in the definitions of these quantities. 

The next step is to solve the renormalization equa- 
tions (2) and (4). To stay within the accuracy of the 
method, in the solution we need to retain the terms 
cx Q:"[ln(A/e)]" while same-order terms with a lower ex- 
ponent of the logarithmic factor should be disregarded. 
This allows for a substantial simplification. We proceed 
by solving Eqs. (2) and (4) at higher energy (far from 
the resonance), where the reflection from the compound 
scatterer is almost perfect. In this case, we approximate 
VL^ni—E)] K, 1. It is possible to sec that in this case the 
renormalization of the tunnel amplitudes t\^2 of each con- 
stituent of our compound scatterer occurs separate from 
each other, dlnti 2/ dine = ai^ ji/2. This renormaliza- 
tion can be incorporated into the energy dependence of 
the effective level widths, TR^L{e) = rfl,i,(e/A)««.^. The 
result for |t(e)P thus reads 
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The above approximation of the integrand in Eq. (2) 
becomes invalid at lower energies, where the transmis- 
sion coefficient may become of the order of unity. The 
energy scale e at which this occurs can be evaluated from 
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Eq. (5), and is given by the solution of equation 2e = 
rL(e) +r/{(e). In the simplest case oi ul = (Xr = a <^ 1, 
it is 2e = (Tl + ^R^i^L + Vr)/2KY. At energies below 
e, the reflection amplitudes in the integrand can be ap- 
proximated as r(e') « 'i'{t). Under this assumption, we 
immediately recover Eq.(3). At |e| < e, its solution yields 



(e - A)2 + f4e)f ^(e) + [F^e") - Tn{^f/A 

(6) 

with 

fLA^) = TL,n{e){\e\lef-^. (7) 

Relations (6) and (7) determine the full crossover func- 
tion for the resonant tunneling between the interacting 
Id electron systems, if e > |A|. 

In the opposite case of a resonance distant from the 
Fermi level, |A| > e, we shall change the approximation 
at e = |A|. The answer is thus given by the equations 
(6), (7) with e being replaced by |A|. 

A typical energy dependence of the transmission prob- 
ability is sketched in the insert of Fig. 1. It combines an 
overall Lorentz-like shape with a sharp dip at the Fermi 
level. Note, that the transmission is not suppressed at 
low energies for a perfectly symmetric barrier with the 
resonant level tuned to coincide with the Fermi level, in 
full agreement with Ref. [7]. 

It may seem that the numerical factors in the defini- 
tion of the crossover energy e and in the condition | A| = e 
of the crossover between low energy cut-offs are chosen 
in a rather arbitrary fashion. Indeed, these two defini- 
tions could contain any other numerical factors of the 
order of 1. The point is that fixing the numerical factors 
with a greater precision would exceed the accuracy of our 
renormalization method. In other words, the energy de- 
pendence of r in all above relations is assumed to be very 
slow. It is this slow dependence that, in the limit a ^ 1, 
gives us the luxury of arbitrary choice of those numerical 
factors. 

To present quantitative conclusions, we discuss the lin- 
ear conductance G{T) in the case of q/v = = Q- 
Within the Landauer formalism, the conductance is given 

by 

G(T) = Gq ^ |t(e)P, (8) 

^ ' ^y_oo 4Tcosh'(e/2T)' ^ ^' ' ^ ' 

where the conductance quantum unit for one fermion 
mode is Gq = e^/27r?i. The results strongly depend on 
the ratio of F/j and F^. We will characterize this ratio 
by the asymmetry parameter (3 = jF^ — Ffl|/(rfl -I- F^) 
which ranges from to 1 and does not depend on energy, 
provided that an = OiL- To emphasize the effect of in- 
teraction, let us recall that in the case of free electrons 



one finds G{T) cx l/T at temperatures T F, A; in the 
limit T ^ 0, the conductance saturates at a finite value, 
which reaches {1 — 0^)Gq if the Fermi level is tuned to the 
resonance (A = 0). Interaction changes this picture no- 
ticeably. Let us start the discussion with the case A = 0. 
At high temperatures, T > e, the conductance can be 
estimated as 

Ga=o(T)/Gq = [7r(l - /32)/4] (T/e)""^ ~ F(T)/T. (9) 

The temperature dependence can be thus ascribed to 
interaction-induced power-law temperature dependence 
of F. Furthermore, the low-temperature behavior dif- 
fers strikingly for symmetric (/3 = 0) and asymmetric 
(/3 ^ 0) resonance. For symmetric resonance, the con- 
ductance saturates at the ideal value of Gq. For /? 7^ 
the conductance reaches at T « e its maximum value, 
which is smaller than (1 — /3^)Gq, and drops to zero with 
the further decrease of temperature, 

Ga=o(T)/Gq = (l//32-l)(r/e~f", T<1. (10) 

The temperature exponents at T ;< e in both cases agree 
with those obtained in Refs. [7,8]. The exponent at /3 ^ 
is the same as for a single tunnel barrier interrupting the 
Id channel. It indicates that at low energies the electrons 
get over the compound scatterer in a single quantum 
transition. The high-temperature exponent arises from 
the separate renormalization of each barrier, which sig- 
nals the "sequential mechanism" of tunneling: the elec- 
tron tunnels across one barrier first and waits a while 
(~ fi,/F(e)) before tunneling across another one. It is im- 
portant to recognize though that no energy relaxation or 
decoherence takes place during this waiting time. This 
is especially clear from our calculation based on the Lan- 
dauer formula: there are no inelastic processes included 
which could provide for the relaxation or decoherence. 

The increase of A leads to a decrease of the conduc- 
tance. For non-interacting electrons, the conductance 
stays at a level of the order of its maximal value, Ga^o 
for A less than F^ -|- T^, which determines the width of 
the resonance in G(A) at T < F^ -|- F/j. At higher tem- 
peratures, the effective resonance width is w ~ T . Let 
us discuss now the temperature dependence w{T) and 
the shape of the resonance G(A) at fixed temperature 
in the presence of interaction. For T S> e, the width 
ui ~ T does not reveal any anomalous exponent. The 
shape of the resonance in this regime is mainly deter- 
mined by the thermal-activated exponential contribution 
G(A) ~ exp(-|A|/T)F(T)/T in Eq. (8). However, at 
large A ::g> w the power-law "cotunneling" tail 

Gtaii(A) = Gq(1 - ff ) {T/if" gVA^ 

replaces that exponential dependence [11]. The crossover 
occurs at A ^ 7'ln(GQ/GA=o) and corresponds to the 
conductance Gcross — G\^q/Gq, this being much smaller 
than Ga=o- 
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At T <C e the apparent width of the non-symmetric 
resonance saturates at w ~ e. The conductance thus 
drops uniformly at any A following the power law (10). 
The symmetric resonance presents an exception. In this 
case, the width shrinks with the decreasing temperature, 

w{T)^iT/er~e, 

and G{T, A) acquires the scaling form, G{T, A) = 
Gq/{1 + [A/w{T)]^}, in agreement with Ref. [7]. 

We further illustrate our results by a numerical eval- 
uation of Eq. (8), see Figs 1-2. For this calculation, we 
choose a = 0.2. By virtue of our approach, the relative 
accuracy of the results is expected to be of the order of a. 
The dependence G(T) is not monotonic, and in the limit 
T ^ the conductance drops to zero at any /3 7^ 0, al- 
though for small f3 this is noticeable only at very low tem- 
peratures ( Fig. 1). The temperature dependence w{T) 
of the width of the resonance G'(A) is shown in the left 
panel of Fig. 2. If /3 ^ 0, this dependence saturates at 
some value w(0) 7^ 0. 
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FIG. 1. Temperature dependence of resonant (A = 0) tun- 
neling conductance. The asymmetry parameter /3 = (top 
curve), 0.2, 0.4, 0.6, and 0.8 (bottom curve). For symmetric 
resonance (/3 = 0), the conductance saturates at T = 0. Inset: 
The typical energy dependence of transmission coefficient. 
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FIG. 2. Left; Half width at half maximum w vs. 
temperature T for the values of asymmetry parameter 
P = 0,0.2,0.4,0.6,0.8 (bottom to top curve). With the 
decreasing temperature, the half width saturates for a 
non-symmetric resonance, and continuously decreases for the 
symmetric one. Right: The conductance dependence on the 
position of the resonant level with respect to the Fermi level, 
G(A), for symmetric (top) and non-symmetric with (3 = 0.5 
(bottom) resonances at three temperatures T/e — 0.04, 0.2, 1. 

The differences and similarities of symmetric and non- 
symmetric resonances are further illustrated in the right 
panels of Fig. 2. The three pairs of line shapes there 
correspond to "high", "medium", and "low" tempera- 



tures, respectively. The two high-temperature curves 
(the smallest values of Ga=o) are hardly distinguish- 
able from each other, and correspond to the resonance 
width w ~ T. Both medium-temperature curves show 
a more narrow resonant peak with increased conductiv- 
ity Ga=o, and are still similar to each other, apart from 
the scale. The real difference becomes visible for the 
low-temperature curves. In the case of non-symmetric 
resonance, the low-temperature curve is just reduced in 
height with no noticeable change of the shape. This is 
in contrast to the symmetric resonance, where the reso- 
nance peak gets taller and thinner. The symmetric res- 
onance does not seem to be very realistic since any ran- 
domness in the barriers and/or in nanowire would cause 
asymmetry. Provided symmetry is achieved, the sym- 
metric resonance can be easily identified by its ideal con- 
ductance Gq. 

To conclude, we have investigated the transmission 
resonances of interacting electrons in Id wires. For a 
weak electron-electron interaction the transmission can 
be considered as an elastic process, which allowed us to 
build a comprehensive theory of the resonances, valid in 
a broad range of temperature and parameters of the res- 
onant level. The temperature dependence of the maxi- 
mum conductance in general is not monotonic, and re- 
veals important differences between symmetric and non- 
symmetric resonances. The obtained quantitative results 
present a comprehensive and consistent picture of the ef- 
fect. It assures us in the qualitative validity of the pic- 
ture at an arbitrary interaction strength. Although we 
are not able to come up with an explicit expression for 
the crossover function G{T) in this case, such function, 
with known high- and low-temperature asymptotes, does 
exist by virtue of the renormalizability. 
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